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We investigate the electronic states of underdoped La2- :E Sr :E Cu04 (LSCO) by using a microscopic 
model, i.e., t-t'-t"-J model, containing vertical charge stripes. The numerically exact diagonalization 
calculation on small clusters shows the consistent explanation of the physical properties in the angle- 
resolved photoemission, neutron magnetic scattering and optical conductivity experiments such as 
the antiphase domain and quasi-one-dimensional charge transport. The pair correlation function of 
the d-channel is suppressed by the stripes. These results demonstrate a crucial role of the stripes in 
LSCO 

PACS numbers: 74.20.Mn, 71.10.Fd, 74.25. Jb, 74.72.Dn 
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Since the discovery of high T c superconductivity, 
La2- 2; Sr 2 ;Cu04 (LSCO) has been extensively studied as 
a typical cuprate superconductor, because it has a sim- 
ple crystal structure with single Cu02 layer and the hole 
density in the Cu02 plane is changeable in a wide range 
from x=0 to 0.35. Understanding of the electronic states 
of LSCO is, therefore, of crucial importance for the study 
of the superconductivity. 

Recently, Ino et al. [0 have performed the angle- 
resolved photoemission spectroscopy (ARPES) experi- 
ment on LSCO and have reported that the spectrum near 
(7r/2,7r/2) along the (0,0)-(7r,7r) direction is very broad 
and weak in underdoped samples. The feature is different 
from the case of underdoped I^S^CaC^Os+a (Bi2212) 
where a sharp peak appears near (tt/2,tt/2) ||. The dif- 
ference between the two typical families may provide a 
clue of high T c superconductivity. Other experimental 
data in LSCO have also revealed anomalous behaviors 
of structural HQ, electronic |p| and magnetic prop- 
erties [Q. In particular, at x=0.12, an incomplete phase 
transition from the low temperature orthorhombic (LTO) 
phase to the low temperature tetragonal (LTT) phase 
(approximately 10% in volume) has been observed ||. 
At the same hole density, the incommensurate antiferro- 
magnetic (AF) long-range order has been reported ||. 
It is interesting that the LTT phase of Nd-doped LSCO 
with £=0.12, Lai.4sNdo.4Sro.i2Cu04, has shown similar 
long-range order accompanied by charge order, which is 
interpreted as charge/spin stripe order that consists of 
vertical charge stripes and antiphase spin domains [fl0|| . 
Therefore, such stripes are expected to play an important 
role in the family of LSCO. From the theoretical side, dis- 
order or fluctuation of the stripe phases has been argued 
as essential physics of high T c superconductors |pd|-[l3|| . A 
phenomenological description of the excitation spectrum 
has also been shown |lj] . 

In this Letter, the electronic states of LSCO are ex- 
amined in terms of the effect of the stripes on various 
excitation spectra. We employ a microscopic model with 
realistic parameters for LSCO (the t-J model with the 



long-range hoppings) under the presence of a potential 
which stabilizes the vertical stripes, simulating the ef- 
fect of the LTT fluctuation considered to be favorable 
for the stripes The numerically exact diagonaliza- 

tion (ED) method is used for small clusters. We find 
that the vertical stripe formation causes the suppression 
of the single-particle excitation A(k, u>), consistent with 
the ARPES data in underdoped LSCO 0. This is a di- 
rect demonstration that the vertical stripes actually exist 
in LSCO. The effects of the stripes on other quantities 
(spin correlation, optical conductivity, and pair correla- 
tion) are investigated and the implication of the results 
is discussed. 

The t-J Hamiltonian with long-range hoppings, termed 
the t-t'-t"-J model, is 

H = J £ Si-Sj-t £ cU, 

E c *W-*" E c t^+H.c. , (i) 

where the summations (i,j) lst , (i>i)2 n d anc ^ (*>j)3rd run 
over first, second and third nearest-neighbor pairs, re- 
spectively. No double occupancy is allowed, and the rest 
of the notation is standard. Recent analysis of ARPES 
data has shown that t' and t" are necessary for under- 
standing not only the dispersion but also the line shape 
of the spectral function Q. In LSCO, we estimated the 
ratio t' /t and t" /t to be —0.12 and 0.08, respectively, by 
fitting the tight-binding (TB) Fermi surface (FS) to the 
experimental one in the overdoped sample [111 on the as- 
sumption that in the overdoped region the FS shape of 
the TB band is the same as that of the t-t'-t"-J model. 
These numbers are slightly smaller than those for Bi2212, 
t'/t=-Q.M and i"/i=0.23 (i=0.35 eV) (2). The contri- 
bution of apex oxygen to the band dispersion may be the 
origin of the difference as discussed by Feiner et al. [^5) . 
By performing an ED calculation of A(k, uj) for a 20-site 
square lattice in the t-t'-t"-J model, we have confirmed 
that the quasiparticle (QP) at k=(7r,0) and (0,7r) exists 
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above the Fermi level in the overdoped region, being con- 
sistent with the observed electronlike FS in the highly 
overdoped sample O . 

In the underdoped region, however, the t-t'-t"-J results 
show sharp peaks along the (0,0)-(7r,7r) direction, being 
inconsistent with the ARPES data (!]]. This implies the 
presence of an additional effect, which may be a stripe 
formation. It is controversial whether the t-J model (also 
t-t'-t"-J model) itself has the stripe- type ground state 
(GS) A possible origin of the appearance of sta- 

ble stripe phase is due to the presence of the long-range 
part of the Coulomb interaction |^0|,|l| and/or the cou- 
pling to lattice distortions. In LSCO, the LTT fluctua- 
tion seems to help the latter mechanism In fact, 
the LTT structure makes Cu-0 bonds along x and y di- 
rections inequivalent, leading to directional distribution 
of carriers through the anisotropy of the Madelung po- 
tential at in-plane oxygen sites and that of the hopping 
amplitude between Cu-O. To model the directional hole 
distribution and the tendency toward the stripe instabil- 
ity as simple as possible, we introduce a configuration- 
dependent "stripe" potential V s . The magnitude of V s 
is assumed to depend on the number of holes nh in each 
column along the y direction of the lattice. In the fol- 
lowing, we use a VT8 x \/l8 cluster with two holes as an 
underdoped system. The maximum number of holes to 
be considered here amounts to three for the final states 
of electron removal A(k, ui). Noting that the 18-site clus- 
ter forms three columns under the periodic boundary 
condition (BC), V s {rih) in each column is assumed that 
V s (0) = V„(l) = 0, V„(2) = -2V, and 7,(3) = -W, be- 
ing V>0. Examples for two configurations are shown in 
Fig. 1(a). V s (rih) behaves like an attractive potential for 
holes independent of distance between holes. When the 
three columns are treated as equivalent, the translational 
symmetry is preserved and thus the size of the Hilbert 
space can be reduced. The symmetry with respect to 
90° rotation is, however, broken as expected. We note 
that the present form of V a (rih) is available only for the 
small cluster with up to three holes (22) . When we treat 
systems with more than three holes, we need to modify 
the form of V s (nh) in order to distinguish, for example, 
between the configurations with a single column filled by 
four holes and with two columns filled by two holes each. 

Let us start on the GS properties of the 18-site cluster 
with two holes as functions of the potential parameter 
V. The total momentum of the GS is always (0,0), and 
a level crossing occurs at V/t=0.6 |23|. For 0<V/t<0.6 
the two holes dominantly form a diagonally bounded pair, 
while for V/t>0.6 both holes are mainly on a single col- 
umn with two-lattice spacing. The latter corresponds to 
the charge stripes expected in LSCO. 

Spin correlation around holes in the GS, C s (a, a') = 
< GS|nl l S l+a -S. 1+Q ,|GS > /N h , is shown in Fig. 1(c). 
Here a and a' denote two sites around a hole following 
the labeling convention shown in Fig. 1(b). is the hole- 
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FIG. 1. (a) Illustration of the stripe potential V s for the 
18-site cluster with periodic BC. The cluster is divided into 
three columns represented by the empty, shaded and slash cir- 
cles. When two (three) holes represented by full circles are on 
a single column, V s = —2V (— 3V). (b) Labeling configura- 
tions used for the calculation of spin correlation around holes 
C 3 (a,a'). The empty circle denotes the position of a hole, (c) 
C s (a, a') of the t-t'-t"-J model with the stripe potential V s on 
a 18-site 2-hole cluster. J/t=0.4, t' /t=-0.12 and t"/t=0.08. 
The discontinuity at V /t=0.6 is due to a level crossing. 



number operator at site i, and Nh is the total number of 
holes. In the stripe regime (V /t>0.6), the spin corre- 
lation for the d and / configurations is AF and strong 
enough to be comparable to the value of —0.33 for the 
Heisenberg model. The spin correlation in the a configu- 
ration is expected to be AF due to the hole motion jl2) . 
Our result shows sizable AF correlation for V/t<0.6. The 
magnitude is reduced when V /t>0.6, since the hole mo- 
tion perpendicular to the stripes is suppressed. However, 
the correlation is still AF for V/i<0.8. When V/t>0.8, 
it becomes weak ferromagnetic (FM). This FM correla- 
tion is, however, due to the periodic BC imposed on the 
cluster that makes the same spin orientation favorable 
in the a configuration. In fact, if we use a 5x4 cluster 
with open BC along the x direction, the spin correlation 
across the stripes is always AF with small magnitude 
(^—0.05 for V/t— 1.0). Such AF correlation is consistent 
with the experimental observation of the antiphase spin 
domains fic|| . 

Figure 2 shows the optical conductivity cm (a;) and 
<j±(lu) (parallel and perpendicular to the stripes, respec- 
tively) on the 18-site cluster with two holes. When 
V/t=l, cr\\(u>) shows large Drude weight indicating one- 
dimensional (ID) behavior, whereas the Drude part in 
<j±(uj) is strongly suppressed. In contrast to ideal ID 
systems |2~i| ], the sizable incoherent weight remains in 
ci|(aj) because the charge carriers in the stripes are not 
completely free from the spin configuration in the spin 
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FIG. 2. Optical conductivity, <j±(uj) and txii(a/), of the 
t-t'-t"-J model with the stripe potential on a 18-site 2-hole 
cluster. J/t=0A, t'/t=-Q.12 and t"/t=0.08. (a) V/t=0. (b) 
and (c) V/t=l. The delta functions (vertical bars) for u/t>0 
are broadened by a Lorentzian with a width of 0.2i (solid 
curves). The Drude weights are shown by solid bars at uj/t—0 
multiplied by 0.4. The negative Drude weight in (b) is due to 
finite-size effect of the cluster. The inset of (c) is the Drude 
weight and integrated conductivity averaged over both the 
directions as functions of V/t. 



domains. The quasi-lD nature along the stripes may ex- 
plain an anomalous x dependence of the integrated <j{uj) 
up to 1.2 eV 0: it increases linearly up to x~0.12 but 
above x~0.12 the value is saturated. At a;~0.12, the hole 
density in each stripe is about 0.5. When more holes en- 
ter into the stripes, the kinetic energy of carriers along 
the stripes decreases since the energy has a maximum at 
the density of 0.5 in the ID chain with strong correlation. 
Thus, the integrated decreases. Holes introduced 

into the spin domains, however, increase the integrated 
value. Therefore, competition between these two effects 
causes the saturated behavior above z>0.12 seen in the 
experiments. 

The inset of Fig. 2(c) shows the total Drude weight D 
and the integrated conductivity K averaged over both the 
directions. We find that the difference K — D, which rep- 
resents the magnitude of the incoherent conductivity, for 
V/t=l is larger than that for V/t=0. The large peak at 
w/t~0.6 in Fig. 2(b), whose energy position is sensitive to 
J, is the origin of the difference. These results are consis- 
tent with the fact that midinfrared absorption at around 
u> ~0.27 eV (~0.7t) in LSCO is enhanced in comparison 
with the spectra of YBa 2 CU 3 6 . 6 and Bi2212 |||. With 
further increasing V/t, the difference K — D decreases. 
This behavior may correspond to the suppression of <t(uj) 
observed in Nd-doped LSCO ||. 

Figure 3 shows A(k,u>) for selected momenta. When 
there is no potential [Fig. 3(a)], the large QP peaks are 
clearly seen below and above the Fermi level at (7r/3,7r/3) 
and (27r/3,27r/3), respectively. For V/t=l [Fig. 3(b)], 
however, there is no distinct QP peak and the spectra 
become more incoherent. At (7r/3,7r/3), a peak is seen 
above the Fermi level. This may be due to the fact that 
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FIG. 3. Single-particle spectral function A(lt, u>) of the 
t-t'-t"-J model with the stripe potential on a 18-site 2-hole 
cluster. J/t=0.4, t'/*=-0.12 and t"/*=0.08. (a), (b), and 
(c) are results for V/t=0, 1, and 2, respectively. The delta 
functions (vertical bars) are broadened by a Lorentzian with 
a width of 0.2t (solid and dashed curves for the electron re- 
moval and addition spectra, respectively). The momentum is 
measured in units of it. (d) The density of states obtained by 
^2 k A(k,ui). The zero energy corresponds to the Fermi level 
defined as the middle point between the first ionization and 
affinity states. 



the stripe potential behaves like an attractive force be- 
tween holes. Here, we note that the reduction of QP 
weight originating from the charge stripes is more or less 
seen at all momenta. In the presence of the potential, 
the coupling between charge carriers in the stripes and 
neighboring spin domains becomes weak (see Fig. 2), and 
the spin correlation across the stripes is also small [see 
Fig. 1(c)]. The implication of the two results is an in- 
crease of the number of configurations that participate 
in the GS. According to the configuration interaction 
picture, such an increase results in the suppression of 
the QP weight and the increase of incoherent structures. 
Therefore, the stripe has a general tendency to make the 
spectrum broad. For V/t=2 [Fig. 3(c)], split-off states 
from incoherent structures are clearly seen at uj/t=0.7 
for (7r/3,7r/3) and (27r/3,0). This is due to the localiza- 
tion of carriers along the direction perpendicular to the 
stripes [^4|. For V/t=l, the tendency to the localization 
may be involved in the low-energy peaks at (27r/3,0), 
but can not be recognized at (7r/3,7r/3). The broadness 
of the spectrum at (7t/3,7t/3), i.e., along the (0,0)-(7r,7r) 
direction, is thus emphasized in contrast to the spectra 
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FIG. 4. Two-hole pair spectral lunction Pi(u) of the 
t-t'-t"-J model with the stripe potential on a 18-site 2-hole 
cluster. J/f=0.4, t'/*=-0.12 and t"/t=0.08. (a) d x 2_ y 2-wave 
symmetry (l=d), and (b) s-wave symmetry (l=s). The delta 
functions are broadened by a Lorentzian with a width of 0.02t. 



for other momenta. As a consequence, our results for 
V/t=l explain the ARPES data We also found 

that the broadness cannot be explained by a diagonal 
stripe potential. Figure 3(d) illustrates the density of 
states (DOS). For large V/t, it shows a gaplike feature 
near the Fermi level, which may be related to the sup- 
pressed DOS at the Fermi level observed in the photoc- 
mission experiment |26[| . 

To see the effect of the stripes on the pairing sym- 
metry of carriers, we calculate the pair spectral func- 
tion Pi(w) for two holes added to the Heisenberg ground 
state |27| , where the hole-pair creation operator is defined 
by N^i 2 J^i e ^'( e ) c i,T c i+c,l ( e ' s are vectors connecting 
nearest-neighbor sites) and A;(e)(=±l) depends on s- 
or c? a .2_ y 2-wave symmetry (l—s or d). For V/t—0, the 
low-energy pair fluctuation is dominated by d-wave sym- 
metry |27| as shown in Fig. 4(a). With increasing V, the 
low-energy weight of Pd{uj) decreases, while P s (w) shows 
enhancement of low-energy fluctuation. This implies that 
the stripes suppress the d-wave pairing accompanied by 
slight enhancement of s-wave channel JlT| . 

In summary, we have investigated the electronic states 
of LSCO by using a microscopic model containing the 
vertical stripes. The tendency toward the formation of 
vertical charge stripes due to the LTT fluctuation was 
modeled by the stripe potential, which was added to the 
t-t'-t"-J model. The realistic values of t' and t" that are 
smaller than those for Bi2212 were used in the exact diag- 
onalization calculation on small clusters. Here, we com- 
ment that the smallness of t' and t" also works favorably 
for the stripe formation [|l9| as compared with Bi2212. 
We found that the model explains the ARPES data with 
suppressed weight along the (0,0)-(7r,7r) direction. Mag- 
netic properties and charge dynamics are consistent with 
experimental data. The consistent results suggest that 
the vertical stripe is an essential ingredient for the ex- 



planation of the physical properties of LSCO. The pair 
correlation function of the d-wave channel is suppressed 
by the stripes, implying the d-wave superconductivity is 
less favorable under the presence of the stripes. This may 
explain why T c in LSCO is so low, as compared with, for 
example, Bi2212. Our results thus clearly demonstrate 
the interesting features of the electronic states caused by 
the stripes, and the uniqueness of LSCO among high T c 
cuprates. 
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